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Abstract. We study the breakdown of the magnetization plateau at the 
magnetization M = Ms/3 [Ms is the saturation magnetization) of the S = 
1/2 anisotropic spin chain with ferromagnetic- ferromagnetic-antiferromagnetic 
interactions. We consider the model with the isotropic ferromagnetic (trimer) 
coupling Jp, and anisotropic antiferromagnetic coupling {Jx = Jy = Jaf ^'Hd 
Jz = AJyij?). For the limit of large 7 = Jp/Jaf, the model is equivalent to 
the 5 = 3/2 XXZ chain with the exchange anisotropy A. There is a phase 
transition between the plateau (small-7) and the no-plateau (large-7) regions. 
This phase transition is of the Berezinskii-Kosterlitz-Thouless type, and we 
determine the phase boundary from the numerical diagonalization data. For 
A = 1, in particular, the phase transition between the plateau and the no-plateau 
regions occurs at the point 7c = 15.4. 



1. Introduction 

Recently, there has been considerable interest in the magnetization processes of one- 
dimensional quantum spin systems. In particular, magnetization plateaues have been 
observed in experimental and theoretical works, for cases such as the 5 = 1/2 
ferromagnetic-ferromagnetic-antiferromagnetic trimerized spin chain lU |^] , the S — 1 
spin chain with bond alternation (and single ion anisotropy) ||^, S p|, pj , the S* = 3/2 
spin chain with single ion anisotropy p , iV-leg spin ladders |9|, |1C|, the = 1/2 



zigzag chain with bond alternation lll^, 12|, and the alternating-spin Heisenberg 
chain with S" = 1/2 and 1 |jl^, This phenomenon stems from the strong 

quantum fluctuation due to the low dimensionality. Extending the Lieb-Schultz- 
Mattis theorem |jl^ , Oshikawa, Yamanaka, and Affleck [Q gave a necessary condition 
for the appearance of the magnetization plateau as a fractional quantized form, 
^(5* — (m)) = integer, where p is the periodicity of the magnetic ground states in the 
thermodynamic limit, S is the magnitude of the spin, and (m) is the magnetization 
per site. Because this is a necessary condition, it depends on the details of the model 
whether the magnetization plateau exists or not, even if their condition is satisfied. In 
fact, there are no magnetization plateau for simple antiferromagnetic Heisenberg spin 
chains (except for zero magnetization (Haldane gap)) 
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In this paper, we study the S — 1/2 trimerized XX Z spin chains in a magnetic 
field h: 

L 

H^J2 hMH3j-2,3j-l{^) + i/3,-1.3,(l)) + JAFi?3,.3, + l(A)] - HM (1) 

J = l 

L 

= Jaf J2 [-7(^^3,-2,3,^1 (1) + i^3,-l,3,(l)) + i?3,,3, + l(A)] " hM 

i=i 

where 

(A) = S^S^, + Spy + AS] SI , (2) 
S" {a = x,y, z) is the spin-1/2 operator at the jth site, Jp — ^Jaf, Jaf > 0, and 

3L 
k=l 

We consider the parameter region — 1 < A < 1. For A = 1, this is a model of 
3CuCl2-2dioxane, for which Ajiro et al measured the magnetization process. Hida 
studied this model using the numerical diagonalization for i = 4, 6, 8 systems. In his 
results, for small 7(7<7c = 2~3) there exists a plateau in the magnetization 
process at M = = Ms/S where Ms = 3L/2 is the saturation magnetization, 

whereas it seems to disappear for 7 > jc- For 3CuCl2-2dioxane, Hida evaluated the 
trimer coupling as 7 = 4.5. Roji and Miyashita jisj calculated the magnetization 
curve by means of quantum Monte Carlo simulation for 7 = 5 and a magnetization 
plateau did not appear. In these studies, however, the phase boundary between the 
plateau and the no-plateau regions is not specified clearly. 




hcl hc2 he he 



Figure 1. Schematic magnetization curve: (a) for 7 < 7c; (b) for 7 > 7c. he is 
the saturation field, (m) is the magnetization per trimer M/L. 

For this model, schematic magnetization curves are shown in figure 1. For the 
weak- ferromagnetic-coupling case (figure 1(a)), the magnetization plateau appears in 
the interval hci < h < hc2, where the magnetization remains 1/3 of the saturation 
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magnetization. In a crude argument, two spins S^j and S^j^i (connected by 
antiferromagnetic bond) tend to form a singlet, and the remaining L spins are aligned 
upward. In the limit 7 — > oo, the low-energy behaviour of the present model effectively 
reduces to that of the S' = 3/2 XX Z spin chain 

L L 

H = J3/2 E [^J^^i + ^Pj+i + ^^l^l+i] - ^ E ^1 (3) 

where J3/2 = Jaf/Q and Sj ~ S'3j_2 + S'3j_i + S'3j, and in this limit it is believed that 
there is no magnetization plateau at M = Afi/3. Thus there should be a transition 
between the plateau and the no-plateau (figure 1(b)) regions at the magnetization 
AI — A/1/3. According to the bosonization approach [0, this transition should be 
of the Berezinskii-Kosterlitz-Thouless (BKT) type g^, |l] , and the width of the 
plateau Ah{M — M1/3) — hc2 — hd near the transition point behaves singularly as 

Ah ^ exp{-C/Vl^) (4) 
where C is a constant. This singular behaviour makes it difficult to study the critical 
properties numerically. 



In our previous paper |22 , we studied an S ^ 1/2 trimerized spin chain: 

L 

H = ^{{1- t)[i/3,-2,3,-l(A) + i/3,-l,3,(A)] + (1 + 2i)i/3,,3, + l(A)} ■ (5) 

i=i 

For this model, we found that in the region surrounded by the three points (A, t) ~ 
(—0.729, 0), (—1, 1), and (—1, —0.5), there is a no-plateau region at the magnetization 
M = Ms/'i. On the basis of conformal field theory and renormalization group 
analysis, we determined this phase boundary successfully. The multicritical point 
(A, i) = (—0.729, 0) is consistent with the numerical Bethe-ansate calculation: (A, i) = 
(—0.729043, 0) [0. In the present study, we apply the same analysis and the numerical 
approach to the model (P. 

The plan of this paper is as follows. In the next section, to study the critical 
properties, we consider an effective continuum model by use of the bosonization 
technique. In section 3, we firstly explain the method used to determine the phase 
boundary between the plateau and the no-plateau regions at Af = Afi/3 numerically. 
Next, we show our numerical results for finite-size systems, and give the boundary 
that we obtained between the plateau and the no-plateau regions. We also check the 
consistency as regards the critical behaviour. The last section is devoted to a summary 
and discussion. 



2. Phase transition 

To see the appearance and the disappearance of the magnetization plateau, we 
introduce the effective-continuum model of the Hamiltonian (|^). In the calculation, 
it is convenient to use the unitary transformation S'fj_i — > —S^^^i, ^tj-i ^ ^^tj-i^ 
S^j_i for all j, and the following parametrization |2j: 

L 

= Jo E [(1 + '5±)(^-2,3,-l + ^-1,3,) + (1 - 2S±H,3, + l] (6) 

i=i 

L 

+ Jo E [(^0 + ^-)(^3.-2^l,-i + + (Ao - 2S.)SI^SI^+,] - hM 
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^f,k — ^j^k ~^ ^j^k C^) 

where 

•^0 — Jaf—^ — 

5± =(7-l)/(27 + l) (8) 
Ao = (-27 + A)/(27 + l) 

Sz - -(7 + A)/(27 + l). 

Using the Jordan- Wigner transformation we firstly map the spin system to a spinless- 
fermion system. Then hnearizing the dispersion relation of the spinless fermions 
around the Fermi point, and with the standard bosonization procedure, we obtain an 
effective-continuum model which describes the low-energy properties of the original 
system. 

The bosonized Hamiltonian of this model was obtained by Okamoto ||^, as 

2' 



H = 



da; 
2^ 



y4,vs I '^cos{V2(t) + 2{kF ~Tr/3)x) (9) 



where vs is the sound velocity of the system, 11 is the momentum density conjugate 
to 0, [(j>{x), n(a;')] = iS{x — x'), and 

TT / M 
fcp = - 1 



2 V Ms, 

is the Fermi point of the spinless-fermion model. We define the lattice spacing as 
1. Due to the oscillation, the second term is important only for M = (i.e., 
fcp = 7r/3), and in the following we consider the system at this magnetization. The 
parameters vs, K, and are related to Jq, S±^z and Aq as 

^ y4,vs = -ttM25^ + S,) (10) 



A=^(l + -^Ao) C^27r(l-^]. (11) 



where 

The dual field 6 conjugate to (p is defined by dxO = ttH, and wc make the identification 
(j) = (j) + \/27r, 9=9 + \/2tt. The spin operator in the continuum picture is given by 

S^{x) = 1 + -^-^(t){x) + I cos(V20 + 2kFX - 7r/3). (12) 
D \/2tt ax 3 

In the Hamiltonian (^) , the first term has a scale invariance and the second term 
violates it. For the free-field case = 0, the scaling dimension of the operator 
•\/2 cos a/20 depends on the Gaussian coupling K as x = K/2. Thus the second term 
of equation (^ is relevant or irrelevant according to whether the renormalized value 
oi K is K < 4 (x < 2) or K > 4 {x > 2). In case where ^ 0, using the notation 
= 4(1 + yo/2) near _fC = 4, we have the following renormalization group equations 

dyo(i) dy0(L) ^T^ ,T\ ^^^^ 

^^ = -y.(L) ^li^--^o(i)l/,(L) (13) 

where L is an infrared cut-off. For yo > yif> > 0, the scaling field y^ is renormalized 
to 0, and yo goes to a finite value; the system is critical and the magnetization 
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plateau does not appear at AI = -A-^i/s- On the other hand, for j/q < 2/0: 2/0 hows 
to infinity and K flows to 0. The phase field is locked to (f> ^ 7r/V2- For the 
spin model (|l|), this means that the second term in equation (|^) makes two spins 
connected by the antiferromagnetic bond form a singlet, resulting in the appearance 
of the magnetization plateau at M — M1/3. The BKT transition occurs at yo — y4>, 
and at the transition point, we obtain 

yo ln(L/Lo) + 1 

where 2/o(> 0) a-nd Lq are some constants. Thus yo{L) flows to very slowly. Near 
the transition point and 2/0 < 2/0: the energy gap — that is, the width of the plateau — 
behaves as equation [ pl| . 

3. Numerical method and results 
3.1. Finite-size behaviour 

In order to study the system numerically, let us consider the finite-size behaviour of the 
model (||). The scaling dimension of the primary field Orn,n — exp(im-\/20 -f- in\/2d) 
for ?/0 = is given by 

Xm,n = y + (15) 

According to the finite-size scaling theory given by Cardy ||2^ , the excitation energy 
of the finite-size system at a critical point is related to the scaling dimension as 

Xm.n{L) = -^{E,nAL) - Eg{L)) (16) 
ITIVS 

where Eg{L) is the ground-state energy of an L-site system with periodic boundary 
conditions (PBC). Near the BKT transition point [K « 4), the excitation energy is 
written as 

^ A£;„,o(i) = ffil„^^ = 2fl + ^^V^ (17) 



2t:vs ' '2 V 2 

-^E,AL) = = - (18) 



^-Kvs 2K{L) 

for integer m and n. In the spinless-fermion language, the integer variable m is 
the difference in number between the left- and the right-moving fermions. The 
variable n is the increment of the fermion number, and relates to the magnetization 
as n = M — Afi/3. Thus considering equation (|l^, we have the logarithmic size 
corrections for finite-size spectrum. 

In the numerical calculation for finite-size systems, we calculate the energy of the 
model (0), fixing the magnetization M = J2'^j with h — 0. If the magnetization 
curve is continuous at the magnetization M = M1/3, the two values h±{Mi/^/ L): 

h+{Ah,^/L) = EiM,/3 + 1,L)~ E{Kh,^,L) 

h-iM,/3/L) = EiM,/3,L)-EiM,/3 - 1,L) (19) 
coincide at L — > 00; also, 

hm /i±(A/i/3/L) = /j(l/2) 
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and this is the magnetic field in the thermodynamic hmit. Here E(M, L) is the lowest 
energy for the L-site system with the magnetization M and h = 0. Then the excitation 
energy in eq. (|l^) can be calculated as 

ASo,n = E{Nh/3 + n,L)~- nh{\/2) - E{My^, L). (20) 

For the critical system, we have Ai?o,n = Ai?o,-n- Averaging these two values we can 
eliminate the magnetic filed as 

Ai?o.±n = i[i?(Mi/3 + +i?(A/i/3 - n,L)] - E{My:,,L). (21) 

In the spinless-fermion language, the magnetic field plays the role of the chemical 
potential 

hiM,/s/L) = [E{M,^3 + 1, i) - ^(Afi/3 - l,L)]/2. 
If the two limits 

lim h+{Mi/j,/L) lim h^{Mi/3/L) 

do not coincide, there exists a magnetization plateau at the magnetization Afi/3, and 
the width of the plateau is given by 

A/i(l/2)= lim [/i+(Afi/3/L)-/i_(Afi/3/L)] = 2 lim ASo,±i(i)- 

To determine the BKT transition point, we apply the method proposed by 
Nomura and Kitazawa p^ , in which a level crossing of some excitation levels is 
used. Under the twisted boundary condition (TBC) S^}^^-^ — —S^'^, S^j^j^i — S^, 
the primary operator Om.n shifts to the operator C'„i+i/2,„ in the quantum spin chain 
[pS] , [2^ , |27[| . For the scaling dimensions of operators \/2 cos((/)/\/2) and \/2 sin(0/\/2), 
we have the following finite-size corrections 

a^±i/2,o(i) = ^ + jyniL) + ^y4,{L) 

^ii/2,o(i) = 1 + - IvAL)- (22) 

The dependence of the coupling comes from the first-order perturbation of the 
second term in equation (^) [ p8| . (Due to the charge neutrality condition in the model 
(||), the operators Oo.n are not corrected in this order.) 

Writing = yQ{l + 1) near the transition point, we have 

xl,/2,o{L) = l-\yo{m + 2t). (23) 
On the other hand, from equation ( p^ ) the scaling dimension of Co, ±2 is given by 

xo,±2{L) = \- iyo(i) = ^A£;o,±2(i). (24) 

From these equations, we see that xq ±2 and 2^5:1/2 cross linearly at the transition 
point (t = 0). 

The energies corresponding to the operators V2cos(0/V2) and svi\{4> / \/2) are 
obtained from two low-lying energies with the twisted boundary conditions as 

^Eli/2fi = ^2:^±i/2,o(i) = E^'''' {Mi/3,L, P = I) - E (A/1/3, L) 



L 

Ai?ii/2,o = ^^ii/2,o(i) = E^'''' {Mi/3,L, P^-l)-E (A/1/3, L) . (25) 
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Here E'^^^ {Mi/^^, L, P ~ ±1) are two low-lying energies under the TBC. The two 
states with these energies are distinguished by the parity of the space inversion 
P : Sj ^ S^L-j+i (in the sine-Gordon model (j) ~^ ~4')- Thus the energy 
differences Aii^o,±2 and AE^^^^ ^ should cross at the transition point. 

3.2. Results 

In the numerical calculation, we consider finite-size systems {L = 4, 6, 8) with the 
PBC and the TBC. Using the above-mentioned degeneracy at the transition point, 
we determine the phase boundary and also check the universality class of the phase 
transition. 

For the PBC, the Hamiltonian is invariant under the translation Sj —>■ and 
the corresponding eigenvalue is the wavenumber q = iirk/L {k = —L/2 + 1, • • • , L/2). 
In the whole region, the lowest energy states with M = Mi/^,Mi/^ ± 2 have the 
wavenumber g = and the parity of the space inversion P = \. 

0.8 
0.6 
0.4 
0.2 

0' 

2 4 6 8 10 

7 




Figure 2. Energy differences A.E'] . (□), Ai?^ , (o), and Ai?o ±2 (+) for 
the 1/ = 6, A = 0.5 system. Tiiese values are defined in equations (Ell) and (ESl) 
in the text. We see a level crossing between AiJ^^^^ g s-iid Ai?o_-|-2. 

Figure 2 shows the behaviour of Ai?o,±2 and Ai?^'^^2 g for L = 6 A = 0.5 systems. 
There is a level crossing between Ai?^j^^2 o ^'^'^ ^^o,±2, as is described in equations 
(p3|), and (p4[). The size dependence of the crossing point is shown in figure 3, from 
which the extrapolated value is estimated as 7c(A = 0.5) = 5.75. 

In figure 4 we show the crossing points for L = 4, 6, 8 on the A — 7 plane. The 
line is the extrapolated value in the polynomials of and this is the boundary 

between the plateau and the no-plateau regions. This line increases with A from 7 = 
at A = —1 to 7 = 15.4 at A = 1. This is explained by the fact that when we decrease 
A from 1, the dimerization effect in the antiferromagnetic bond is weakened. 

As a consistency check with the critical theory, we calculate the following averaged 
scaling dimension: 

^±i/2,o(^) + 3a^±i/2,o(-^) , s 

4 ■ 

Taking this average, we can eliminate the leading logarithmic correction ( p^ ) for finite- 
size systems at the critical point i = (see equation (p3|)). To calculate the scaling 



Breakdown of the magnetization plateau of a spin chain 
6.3 




_i I I I I 

0.01 0.02 0.03 0.04 0.05 0.06 
-2 



1/L' 



Figure 3. The size dependence of the crossing point 7^ for A = 0.5. 
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Figure 4. Crossing points for L = 4 (+), L = 6 (o), and L = 8 (•). The solid 
line shows the extrapolated values. 



dimension, we need the sound velocity vs- This can be calculated using the lowest 
energy with M = A/1/3 ^-nd the wavenumber q = 2tt/L (corresponding to the U(l) 
current): 

E{Ms/S, L, q = 27r/L) - E{Ms/S, L) 



Vs = lim 

L— >oo 



2ti/L 



(27) 



Figure 5 shows the extrapolated value of the averaged scaling dimension (g6|) and the 
bare value of the i = 8 system on the transition line. The averaged scaling dimension 
is very close to the expected value x = 1/2, whereas the bare values o(^) ^^"^ 

^±1/2 o(^) from X — 1/2 due to the logarithmic size correction. 

To confirm the universality class of the phase transition for the M = 
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systems, we also calculate the conformal anomaly number c from [ p9[ pO[ 

^^^^=e(AW3)-^ + .... (28) 

For the BKT critical point, the conformal anomaly number is c = 1. We see that the 
value c is close to 1 within the error of a few percent {e.g. c = 1.01 for A = 1). Thus 
we can conclude that the transition between the plateau and the no-plateau region at 
M = Ml/3 is of the BKT type. 
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Figure 5. Extrapolated values of the averaged scaling dimension (aJ^^^j o ~^ 
■^^±1/2 0^/^ on the transition line (•). We also show the bare values of 
^±l/2,o(-^ = 6. °) and x'^-^^^fi^^ = 6. °)- 



4. Summary and discussion 

We studied the plateau-no-plateau transition of the model (|l|) for the magnetization 
M = Ms/3. For small 7 there is a plateau in the magnetization process, while for 
large 7 there is no plateau (figure 1). By use of the effective-continuum model and the 
renormalization group method, we showed that this quantum phase transition is of 
the Berezinskii-Kosterlitz-Thouless type. On the basis of conformal field theory and 
renormalization group analysis, we argued for the finite-size spectrum and determined 
the phase boundary (figure 4) using a level crossing of special excitations (figure 2). 
This boundary runs from (A, 7) = (—1,0) to (A, 7) = (1,15.4). We also checked the 
consistency of the critical theory and the numerical calculation, and concluded that 
the phase transition is of the BKT type. 

The extrapolated critical point for A = 1 is 7 = 15.4, which is somewhat larger 
than the value from previous experimental and numerical studies [|l|, |l^, |8|. One 
possible explanation for this discrepancy is as follows. The width of the plateau 
near the transition point is given by equation (^), and the coefficient in exponential, 
C, is large for A = 1. Thus the narrow-plateau- width region (that is, the small- 
gap region) is widely ranged. This means that it is difficult to estimate the plateau 
width Ah directly from the numerical calculation and also at finite temperatures. 
According to Kosterlitz [^, the coefficient C is described as C cx and from 
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figure 5, we see that yo{= y,p) is small at tlie transition point for A = 1. For 
the S = 3/2 antiferromagnetic Heisenberg chain corresponding to the large-7 limit, 
the Gaussian coupling K is calculated as if « 4.4 (the compactification radius is 
R = (2ttK)~^/^ = 0.95/V87r) in reference [Q, which is slightly larger than the critical 
value K ^ 4. This fact supports a large value of 7c for A = 1. 
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